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ABSTRACT 



We report the results of a spectroscopic and polarimetric study of the massive, 
0-i \ hydrogen-rich WN6h stars R144 (HD 38282 = BAT99-118 = Brey 89) and R145 (HDE 

269928 = BAT99-119 = Brey 90) i ntheLMC. Both star s have been suspecte d to be bi- 
naries by previous studies (R144: ISchnurr et aUfeOOSbl ; R145: lMoffa3ll989f) . We have 
combined radial-velocity (RV) data from these two studies with previously unpub- 
lished polarimetric data. For R145, we were able to establish, for the first time, an 
orbital period of 158.8 days, along with the full set of orbital parameters, including the 
inclination angle i, which was found to be i = 38° ±9°. By app lying a modified version 



■ of the shift-and-add method developed bv lDemers et alj ( 20021 ). we were able to isolate 



the spectral signature of the very faint-line companion star. With the RV amplitudes 
of both components in R145, we were thus able to estimate their absolute masses. We 
find minimum masses MwR,sm 3 i = (116 ± 33) M Q and Mosin 3 i = (48 ± 20) M Q for 
the WR and the O component, respectively. Thus, if the low inclination angle were 
correct, resulting absolute masses of the components would be at least 300 and 125 
M Q , respectively. However, such high masses are not supported by brightness consid- 
erations when R145 is compared to systems with known, very high masses such as 
NGC3603-A1 or WR20a. An inclination angle close to 90° would remedy the situa- 
tion, but is excluded by the currently available data. More and better data are thus 
required to firmly establish the nature of this puzzling, yet potentially very massive 
and important system. As to R144, however, the combined data sets are not sufficient 
to find any periodicity. 



Key words: binaries: spectroscopic - stars: early-type - stars: fundamental param- 
eters - stars: Wolf-Rayet 



1 INTRODUCTION 

Fundamental questions of stellar astrophysics suffer from 
a lack of truly empirical evidence when it comes to stars 
with the highest masses both on and evolved off the 
main sequence (MS). Recent models maintain that in the 
early Universe, the initial-mass function (IMF) was top- 
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heavy and that primordial (Population III) stars were ex- 
tremely massive, w i th ma s ses from ^100 to ^1000 M^ 
(lOstriker fc Gnedinl Il996l : iNakamura fc Umemural 1200 ll : 
ISchaererll2002l ). In contrast to that, recent studies have put 
forward both theoretical and observational arguments that 
under prese nt-day conditions, an IMF cut-off occ urs around 
~ 150 M Q l|Weidner fc Kroupall20o4 lFigerll2005r i . However, 
even in the Local Group very massive stars remain poorly 
characterized, and it so far remains uncertain at what mass 
an upper cut-off really occurs, if at all. 

The least model-dependent and thus most reliable way 
to directly measure stellar masses is by Keplerian orbits 
in binary systems. In the past years, studies of very mas- 
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sive binaries have yielded the somewhat surprising results 
that the most massive stars known are to be found not 
among O-type stars, but among a very luminous subtype 
of Wolf-Rayet (WR) stars, the so-called WN5-7h (or ha) 
stars. These hydrogen-rich WN stars are young, unevolved 
objects rather than "classical" WR stars, which are usu- 
ally identified with e v olved, helium-burning, ma ssive stars 
l|de Koter et al J 1 19971 ; ICrowther fc Dessardll998h . WN5-7h 
stars mimic classical, i.e. evolved, core helium burning WR 
stars because their extreme luminosity, a result of their very 
high masses, drives a dense and fast wind which gives rise 
to WR-like emission lines. Due to the rarity of very massive 
stars in general and of very massive binaries in particular, 
only few massive systems have been "weighed" so far. In 
order of increasing mass, these are the Galactic Wolf-Rayet 
stars WR22 (72 for the WN7h compone nt (|Rauw et all 
1 19961 ). although ISchweickhardt et al.1 (1 19991 ) derive only 55 
M ) , WR20a (83 and 82 M for both W N6ha compo- 
nents; ^^m^^MoM^m^s^^^Q^), WR21a (min- 
imum mass M sin 3 i = 87 Mq for the 03H7WN6 primary; 
<|Niemela et al.ll2008l )). NG C3603-A1 (116 and 8 9 M for 
both WN6ha components; ISchnurr et alj [2008al V Clearly, 
these very high masses impressively underline the extraor- 
dinary nature of WN5-7h stars. 

In a recent, i ntense survey of 41 of the 47 known WNL 
stars in the LMC. ISchnurr et al.1 l|2008bl . hereafter S08) used 
2m-class telescopes to obtain repeated, intermediate-quality 
(R=1000, S/N~80) spectra in order to assess the binary sta- 
tus of the 41 targets from radial-velocity (RV) variations. 
S08 identified four new binaries containing WN5-7h stars, 
bringing the total number of known WNLh binaries in the 
LMC to nine. For one of the previously known binaries, R145 
(= HDE 269928 = BAT99-11 9 = Br e y 90), the preliminary 
25.4-day period reported by iMoffatl jl989l . hereafter M89) 
could not be confirmed, but neither could a coherent pe- 
riod be established from the RV data of this study alone. 
Therefore, we combined S08's RVs with those published by 
M89, and with previously unpublished polarimetry. It was 
only after the combination that we were able to find the 
true orbital period of this system as well as other orbital 
parameters. 

From both its RV and EW variations and its high X- 
ray luminosity, S08 identified another WN6h star, R144 (HD 
38282 = BAT99-118 = Brey 89), to be a binary candidate, 
although M89 had identified this star as probably single. 
In terms of spectral type, R144 is almost a perfect clone 
of R145; however, at the ~same bolometric correction and 
after allowing for differential interstell ar extinction, R144 is 
^0.5 m ag brighter than R145. Indeed. [Crowther fc Dessarj 
(1998), using atmosphere models without iron-line blanket- 
ing, derived a spectroscopic mass in excess of 100 Mq , and 
a luminosity of logL/Z/Q = 6.34 for R144, which makes 
this star one of the most luminous WR stars known. New, 
updated models most likely would yield an even higher 
mass and luminosity, making R144 the most luminous main- 
sequence (MS) object in the Local Group. Since we also have 
unpublished polarimetry for R144, we have revisited this ob- 
ject as well, using the spectroscopic data obtained by S08. 

In the present paper, we will report the results of this 
study. In Section [2j we will briefly describe the observations 
and the data reduction. In Section O we will describe the 



analysis of our data and present the results. The paper is 
then summarized in Section [4] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectroscopy 

The spectroscopic observations are described in detail in 
S08. We briefly recapitulate here that long-slit spectro- 
graphs attached to ground-based, southern, 2m-class tele- 
scopes were used. Data were obtained in three observing 
seasons between 2001 and 2003 to maximize the time cover- 
age, and were carried out during 13 runs at 6 different tele- 
scopes. The spectral coverage varied from one instrument to 
another, but all spectrographs covered the region from 4000 
to 5000 A, thereby encompassing the strategic emission line 
He n A4686 

Data were reduced in the standard manner using 
NOAO-IRAF0, and corrected for systematic shifts among 
different observatories (see S08 for more details). The 
achieved linear dispersion varied from 0.65 A/pixel to 1.64 
A/pixel, but all data were uniformly rebinned to 1.65 
A/pixel, thereby yielding a conservative 3-pixel resolving 
power of R ~ 1000. The achieved signal-to-noise (S/N) ratio 
was ~ 120 per resolution element for both stars, as measured 
in the continuum region between 5050 and 5350 A. 

Additional RV data for R144 and R145 were taken from 
M89, to which we refer for more details on the data acqui- 
sition and reduction. 



2.2 Polarimetry 

Linear polarimetry in white light was obtained during a to- 
tal of six observing runs between October 1988 and May 
1990 at the 2.2m telescope of the Complejo Astronomico El 
Leonc ito (CASLEO) near Sa n Juan, Argentina, with VAT- 
POL (|Magalhaes et al.l Il984h . and at th e ESO/MPG-2.2m 
telescope at La Silla, Chile, with PISCO l|Stahl et al.lll986l ). 
RCA GaAs phototubes were employed in both cases. Expo- 
sure times were typically 15 minutes per data point. Appro- 
priate standard stars were taken to calibrate the polarization 
angle and the zero level, and exposures of the adjacent sky 
were obtained to correct for the background count rates. 

The polarimetric data w ere calibrated the usual way 
(e.g. IVillar-Sbam et al.l [2006) . Since Thompson scattering 
is wavelength-independent (i.e. grey), and the differences 
between the sensitivity curves of the respective detectors 
are small, no correction for the different passbands was 
deemed necessary. However, a small instrumental offset be- 
tween PISCO and VATPOL was found in polarization; we 
therefore have added AQ = 0.10% and AU = 0.10% to 
the PISCO data to match the VATPOL data, on average. 
The final polarimetric data for R144 and R145 are listed in 
Tables [T] and respectively. 

Statistics reveal that R144 displays a small but signif- 
icant scatter (standard deviation) about their mean values 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Linear polarimetry for R144, obtained with VATPOL 
at CASLEO over a two-week period in January 1990. 



HJD 


Q 


°"<3 


U 




Instrument 


2447907.783 


-0.292 


0.022 


0.175 


0.022 


VATPOL 


2447908.819 


-0.213 


0.031 


0.112 


0.031 


VATPOL 


2447909.813 


-0.103 


0.028 


0.099 


0.028 


VATPOL 


2447910.681 


-0.136 


0.020 


0.206 


0.020 


VATPOL 


2447911.754 


-0.131 


0.026 


0.204 


0.026 


VATPOL 


2447912.754 


-0.084 


0.025 


0.160 


0.025 


VATPOL 


2447913.566 


-0.162 


0.020 


0.188 


0.020 


VATPOL 


2447914.628 


-0.134 


0.035 


0.213 


0.035 


VATPOL 


2447915.762 


-0.138 


0.021 


0.220 


0.021 


VATPOL 


2447916.758 


-0.103 


0.028 


0.232 


0.028 


VATPOL 


2447917.803 


-0.194 


0.030 


0.218 


0.030 


VATPOL 


2447919.763 


-0.179 


0.031 


0.308 


0.031 


VATPOL 


2447920.753 


-0.173 


0.025 


0.277 


0.025 


VATPOL 


2447921.699 


-0.180 


0.025 


0.279 


0.025 


VATPOL 


2447923.789 


-0.170 


0.025 


0.310 


0.025 


VATPOL 
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Figure 1. Polarimetric data for R144 plotted versus Julian date. 

in both Stokes parameters, 0.052% and 0.063% in Q and U, 
respectively; by comparison, the quoted error per data point 
is only ~half this. However, in Figure [JJ one can clearly see 
that over the observation period of ~15 days, the polariza- 
tion in Stokes-(7 increases constantly by ~0.2% (i.e. ~seven 
times the quoted error per data point), while Stokes-Q also 
shows systematic variability, which is a strong indication 
that some coherent process is involved, e.g. a (longer-period) 
binary motion. 

R145, on the other hand, is much more variable; the 
scatter is 0.26% and 0.20% about their mean values, respec- 
tively, in Q and U , with well defined and coherent variability 
episodes (Figure [2). 

3 DATA ANALYSIS AND RESULTS 
3.1 Spectroscopic Data 

We here only briefly describe the methods employed to ex- 
tract radial velocities (RVs) from the spectra (for details 
we refer the reader to S08 and references therein): Cross- 
correlation and emission-line fitting routines programmed 



Table 2. Linear polarimetry for R145, obtained with PISCO at 
ESO and VATPOL at CASLEO, between October 1988 and May 
1990. 



HJD Q ctq U ay Instrument 



using the ESO-MIDAS tasks XCORR/IMA and FIT/IMA, 
respectively, with an iterative approach for cross-correlation. 

Since the spectra of both R144 and R145 are largely 
dominated by the strong He n A4686 emission line (see Fig- 
ure [3j, we initially confined cross-correlation to this line to 
derive RVs for R145, but then applied the method to the 
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Figure 2. Polarimctric data for R145 plotted versus Julian date. 
Note the enormous amplitude of rarely seen so large in 

other WR+O systems. 



Hell A5412 line as well. Furthermore, RVs from the faint, 
but very narrow Niv A4058 emission were obtained by fit- 
ting Gauss profiles to the line. 

Unfortunately, due to low S/N and an unreliable wave- 
length calibration at the blue end of our spectra (the 
comparison-arc lamp has practically no useful lines in this 
region; see S08 for more details), the Niv A4058 line yielded 
much worse RV scatter than the two Hell lines. Thus, the 
main part of the analysis was carried out with the two, much 
stronger He II lines 

Fitting a Gauss function to the Hell A4686 emission 
yielded a slightly larger scatter than cross-correlation, but 
was used to measure absolute, systemic velocities. These ab- 
solute velocities were then used to correct the relative RVs 
obtained by cross-correlation, and to combine our data with 
those of M89. Initially, a possible zero-point shift between 
S08 and M89 was not corrected for (but see below). Unfortu- 
nately, M89's spectra did not cover HellA5412, so that only 
RV data from Hen A4686 and Niv A4058 could be com- 
bined. The combined RV data for Hell A4686 are shown in 
Figure 2] 
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Figure 3. Mean spectra of R144 and R145. For clarity, the upper 
spectrum has been shifted by 2 flux units. 



3.2.1 R144 

On close inspection in Figure 2] R144 does not show any 
clearly repeating pattern in S08's RV data, which is the rea- 
son why S08 were not able to establish a period in the first 
place, and there certainly is no pattern at all visible in the 
M89 data set. We therefore have to conclude that better 
data are required for R144, if any periodicity is to be found. 
Long-term monitoring is being undertaken. 



3.2 Search for Periodicities in the RV Curves 

As reported by S08, a comprehensive period analysis us- 
ing different methods of their RV data in the range from 
1 to 200 days yielded no coherent period for R145. We re- 
peated this analysis on the combined S08+M89 data set, also 
without success. Additionally, we tried t o apply the phase - 
dispersion minimization (PDM) method (|Stellingw crf 1978), 
which is very reliable because it does not assume a specific 
wave form (e.g. sinusoids) and which is therefore very use- 
ful for strongly eccentric binaries with highly distorted RV 
curves. Unfortunately however, the data are too dispersed 
for periods much longer than 100 days, so the PDM suffers 
from many severely depleted or downright empty phase bins. 
Thus, a different approach was used. 



3.2.2 R145 

Inspection of R145's RV curves in Figure [3] reveals a repeat- 
ing pattern: In both M89's and S08's last observing season, 
the RV curve drops from its highest to its lowest value within 
a very short interval of time. This event can consistently be 
interpreted as the periastron passage (PP) of a highly ec- 
centric, long-period binary. But what is its period? We are 
lucky, because it seems that we have three such PP events 
recorded: two in the M89 data set (called PP1 and PP2), 
and one in the S08 data set (PP3). We further distinguish 
the highest and the lowest RVs during the respective event; 
thus, we find that the combined RV data set covers the fol- 
lowing events: PP1 high, PP2 high and low, PP3 high and 
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Figure 4. RV data obtained from Hell A4686 of R144 (left) and R145 (right). Shown are the three observing seasons of M89 (upper 
panels) and S08 (lower panels), respectively. 



Table 3. Measured dates and RVs for the presumed periastron- 
passage (PP) events PP1 through PP3 in the data sets of M89 
and S08, respectively. "High" and "low" refer to maximum and 
minimum RVs measured just before and after the supposed PP, 
respectively. See text for more details. 



Data set 


Event 


MJD a 


RV 6 








[kms- 1 ] 


M89 


PP1 high 


44261.4 


423 ±15 


M89 


PP2 high 


44584.4 


424 ±15 


M89 


PP2 low 


44599.5 


276 ±15 


S08 


PP3 high 


53001.8 


432 ±15 


S08 


PP3 low 


53012.4 


252 ±15 



a JD-2,400,000.5 
'errors are estimated 



low. Julian Dates and RVs for each event are given in Table 

El 

It seems that the duration of the PP2 event is somewhat 
longer than that of the PP3 event, ~15 versus ~10 days. 
However, the observed RV amplitude during PP2 is smaller 
than that observed during PP3. It is thus likely that the 
true minimum RV was missed during the observations of the 



PP2 "low" event, thus the full RV peak-to-valley swing of 
the WR star (i.e., 2K) might be closer to the value obtained 
from PP3, ARV(PP3) = 180 kms" 1 . 

Measuring the elapsed time between corresponding 
events yields the following time intervals, each of which will 
be an integer multiple of the period, assumed, of course, 
there is one: 

AT 31 , M = MJD(PP3 hi - PPlhi) = (8740.4 ± 6)days 
AT 32 ,M = MJD(PP3 hi - PP2 hi ) = (8417.4 ± 6)days 
AT 32 ,io = MJD(PP3io - PP2i ) = (8412.9 ± 6)days 
AT 2 i, hi = MJD(PP2 M - PPl h i) = (323.0 ± 6)days, 

where the overall, quadratic error has been obtained from 
\/ A 2 + 4 2 ~ 6. From AT2i,hi follows that the longest possible 
period is (323 ± 6) days; otherwise, two high events would 
occur during one orbit. The true period is thus an integer 
fraction of this period, i.e. (323 ± 6)/n days, with n = 1..5, 
since any period shorter than ~60 days is ruled out by S08's 
observations; it follows that there are five groups of periods, 
323±6, 161. 5±3, 107.67±2, 80.75±1.5, and 64.6±1.2 days. 
Within each group, the longer time intervals are used to 
refine the results and to obtain more precise periods. Periods 
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Table 4. The five groups of possible orbital periods of R145 as 
computed from the condition that the true period has to be an 
integer fraction n of the respective time interval AT2i i hi, AT32J01 
AT32 hii an d AT 31 hi . For each group, the typical error on the 
period is given. Periods that are consistent with each other within 
these errors, have been averaged to yield P mean and the mean 
error. 





Possible periods 


(in days) 




P21M 


P32,lo 


P32,hi 


Psi.hi 


fmcan 


323.00 


323.573 


323.746 


323.719 


323.679 


±6.00 


±0.231 


±0.231 


±0.222 


±0.228 


161.50 


161.787 


161.874 


161.859 


161.840 


±3.00 


158.734 


158.819 


158.916 


158.823 




±0.115 


±0.115 


±0.111 


±0.114 


107.67 


106.492 


106.549 


106.590 


106.544 


±2.00 


107.858 


107.915 


107.906 


107.906 




109.258 


109.317 


109.255 


109.277 




±0.078 


±0.078 


±0.075 


±0.077 


80.75 


79.367 


79.409 


79.458 


79.411 


±1.50 


80.123 


80.166 


80.187 


80.159 




80.893 


80.937 


80.930 


80.902 




81.679 


81.722 


81.686 


81.696 




±0.058 


±0.058 


±0.056 


±0.057 


64.60 


63.734 


63.763 


63.799 


63.765 


±1.20 


64.221 


64.255 


64.268 


64.248 




64.715 


64.749 


64.744 


64.736 




65.216 


65.251 


65.227 


65.231 




65.726 


65.761 


65.717 


65.735 




±0.047 


±0.047 


±0.043 


±0.046 



which, within their errors, are consistent with each other, 
are averaged, and a quadratic mean error is computed. The 
resulting periods are listed in Table [3] 

Since the data are too sparse, an autom ated phase- 
dispersion minimum method l|Stellingwerflll978l ) is not pos- 
sible; we thus folded the data into the respective phases and 
chose the most coherent-looking curve by eye. This method 
is of course limited; for instance, there are no noticeable 
differences in the RV curves of the 323-day and 161.5-day 
periods within a given group. Hence, the selected RV curves 
shown in Figure [5] are representative. However, the groups 
of periods around 107.67, 80.75, and 64.60 days can readily 
be ruled out as they yield incoherent RV curves. 

The longest possible (and unique) period, 323.679 days, 
results in a very well-defined RV curve and a very eccentric 
orbit; however, it does so by construction and, more im- 
portantly, shows a data gap of almost exactly half a phase, 
which is suspicious. The periods around 161.5 days also yield 
coherent RV curves, but with a much better phase coverage, 
and a slightly less eccentric orbit, as can already be seen by 
eye. 



3.3 Polarimetric Data 

With these results from the spectroscopy, we proceeded to 
analyze the polarimetric data. To derive the o rbital param- 
eters, we used the elliptical or bit model by Brown et al.l 
(1982) with the correction by ISimmons fc Bovld i| 19841 ). 
modified for an extended source of scatterers (cf. 
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Figure 5. Combined RV data of R145 obtained from He II A4686 
(M89: filled; S08: empty symbols) and folded into the correspond- 
ing phases of five different, possible periods based Oil Pmean 1H 
Table H 



iRobert et ailll992l ). i.e. assuming an ensemble of optically 
thin, free electrons that follow the wind density around 
one of the stars s pherically symmetrically (for details, see 
iMoffat et al]|l998h . 

For R144, the fit did not converge, and since there is 
no orbital solution from the spectroscopy, we were unable to 
constrain the fit. This does not entirely come as a surprise, 
given the small number of data points (N = 15) spread 
over only ~2 weeks, and the rather small level of variability. 
Thus, we abandon the study of R144 based on the current 
data. 

For R145, however, the fit of the polarization data alone 
yielded a coherent solution with a period of P = (159 ± 1) 
days, which is consistent with the lower boundary of P2i,M 
(see above), but has consequences that will be discussed 
below. Encouraged by this result, we combined the spec- 
troscopic (i.e., M89±S08) and polarimetric data sets, and 
forced a simultaneous fit. For cross- verification, we also used 
the RV data that were obtained from the slightly weaker, 
but less perturbed, Hen A5411 emission line, with similar 
results. 

Note that the fits to the respective data sets do not 
have all orbital parameters in common. For instance, spec- 
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troscopy cannot yield the inclination angle i and the orien- 
tation of the line of nodes both of which are calculated 
from the polarimetry. Thus, only the following parameters 
were obtained from the fit of the combined RV and po- 
larimetric data: Interstellar polarization Qo and Uo; (see 
above); the orbital inclination angle i; the orbital longitude 
of the centroid of the scattering region at periastron passage 
X p (linked to the argument of the periastron by the relation 
\ p — uj — 90°); the orbital period P; the RV amplitude -Kwr 
of the WR component; the systemic velocity Vb; the system- 
atic RV shift between the M89 and S08 data sets AV; the 
orbital eccentricity e; the time of periastron passage To; the 
total electron scattering optical depth r (assuming Thomp- 
son scattering to be the only source of opacity) ; the exponent 
7 in the power law describing how the electron density falls 
off from the WR stars (e.g., 7 = 2 for an inverse square law 
around a scattering point source); the possible systematic 
shift between the two RV data sets ARV = S08 - M89. 

In order to facilitate convergence, the fit was carried out 
for three parameter groups which are, to first order, inde- 
pendent of each other: 1. (K, Vb, ARV); 2. (Qo, Uo, r, 7); 3. 
(i, fl, X p , e). One group was fitted while the two others were 
kept constant. For each permutation x 2 was computed, and 
the procedure was repeated until the fit yielded a constant 
X 2 for all permutations. In order to estimate the errors on 
individual fit parameters, we explored the sensitivity of the 
fit (i.e., its % 2 ) by varying one parameter around its best-fit 
value (obtained from the overall fit) while keeping all other 
parameters constant. Those variations of each individual pa- 
rameter which resulted in a 5% deterioration of the x 2 com- 
pared to its minimum value from the overall fit, were used 
as error levels of the respective parameter. Results of the 
fits for both Hen A4686 and Hen A5411 are shown in Table 
[5] The resulting orbital solutions together with the polari- 
metric data folded into the corresponding phase is shown in 
Figure [6] 

The difference in the systemic velocities obtained from 
Hen A4686 and A5411 amounts to ~55 kms -1 . It is well 
known in the literature that systemic velocities obtained 
from emission lines of WR stars display a systematic red- 
shift with respect to the true systemic velocity (e.g. of 
the host galaxy of the WR star, see e.g. S08). This phe- 
nomenon ca n be explained by radiative-transfer effects (cf. 
lHillierlll989h . However, it is noteworthy that different lines 
of the same ionic species will display significantly different 
systemic velocities. Cross-verification with R144 yielded the 
very same effect with roughly the same shift (~50 kms - ) 
between the systemic velocities obtained from the two He n 
lines, and shows that the effect is indeed real. 

However, in WR+O binaries, Hell emission lines are 
prone to suffer from profile variations due to excess emission 
arising in the wind-wind collision (WWC) region; therefore 
they might not yield the correct, but a more or less distorted 
RV curve (and with this, wrong orbital parameters). Indeed 
R145 displays strong excess emissions, and because R145 is 

2 We here follow the definition of lHarries &: How arth (1996|), i.e. 
that Q is the angle of the ascending node measured from north 
through east with the constraint that f2 < 180° , i.e. an ambiguity 
of 180° exists in the determination of f2 since it is impossible from 
polarimetry to discern the ascending from the descending node. 
This is done using the spectroscopic orbit of the binary. 



Table 5. Orbital parameters from the combined fit of both the 
spectroscopic and the polarimctric data sets. 



Parameter Hen A4686 Hen A5411 



T [%1 


0.048 


± 


0.006 


0.049 


± 0.006 


Qo 1%) 


-2.17 


± 


0.04 


-2.18 


± 0.04 


Uo [%] 


0.67 


± 


0.04 


0.68 


± 0.04 


n [°] 


-40 


± 


7 


-42 


± 7 


i [°] 


37 


± 


7 


39 


± 6 


lj = \ p + 90 [°] 


85.0 


± 


6.7 


83.4 


± 6.6 


P [days] 


158.8 


± 


0.1 


158.8 


± 0.1 


K W R [kms" 1 ] 


81 


± 


21 


93 


± 21 


ARV a [kms" 1 ] 


10 


± 


29 




n/a 


e 


0.70 


± 


0.02 


0.69 


± 0.02 


Vb [kms" 1 ] 


328 


± 


14 


382 


± 13 


T [JD-2,450,000.5] 


3007.763 


± 


0.25 


3007.763 


± 0.25 


7 


1.70 


± 


0.15 


1.75 


± 0.14 



a ARV=S08-M89 



a very eccentric system, the strength of the excess emission 
varies considerably over the orbital phase. From simple con- 
siderations, we expect (and find; see Section [3.61 for a more 
detailed discussion) the excess emission to be strongest near 
periastron. This, however, means that the RVs derived from 
line profiles of both Hell A4686 and A5411 are most affected 
at an orbital phase where the RV data points have the most 
importance for the orbital fit, namely when the RV of the 
WN component changes quickly from maximum redshifted 
to maximum blueshifted velocities, around periastron pas- 
sage. 

To somewhat reduce the effect of such line-profile vari- 
ations over the orbital cycle, we have employed the iter- 
ative cross-correlat ion method described in S08. In fact, 
iFoellmi et ail (2003) have shown that using bisectors to mea- 
sure the line center yields perfectly similar results to those 
obtained by iterative cross-correlation (with bisectors yield- 
ing a slightly larger scatter). Thus, the iterative form of 
cross-correlation is relatively robust against the influence 
even of strong and variable WWC excess emissions. 

Furthermore, we have investigated how the emission- 
line strength changes over the orbital cycle due to the vari- 
able strength of the excess emissions. We measured the 
equivalent widths (EWs) of both Hell A4686 and A5411, 
normalized them by the mean value around apastron (phase 
interval [0.4;0.6]) when the WWC excess emissions suppos- 
edly are weakest, and folded these relative EWs into the 
orbital phase (for P — 158.8 days; see Figure 0. Despite 
the fact that Hell A4686 suffers more from the increased 
WWC excess emission at periastron (almost twice compared 
to Hell A5411), the orbital solutions the respective lines are 
fully consistent with each other, the only difference being 
the redshifted systemic velocity of Hell A5411, as discussed 
above. 

As an additional test, we have used RVs obtained 
from Gauss-fitting the weak and narrow, highly-ionized N IV 
A4058 emission line. As before, our results were combined 
with those published by M89 using the systematic shift be- 
tween the two data sets obtained from the orbital fitting 
of Hen A4686 (ARV of 10 kms" 1 ; see Table EJ. Due to- 
fts flu orescent origin from ultra-violet transitions (cf. iHillierl 
1988), the Niv A4058 line is commonly regarded as very ro- 
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Figure 6. Orbital solution for R145 obtained from the simultaneous fit of the the Stokes Q and U, and RV data sets. For the Hell 
A4686 emission line (bottom panel), S08 data are shown in filled symbols, and M89 data arc shown in empty symbols. The systematic 
shift between the two RV data sets, ARV, has been applied. For Hen A5411, only SOS data could be used for the fit, which is why there 
are fewer data points. 



bust against WWC excess emissions, and usually yields RVs 
which follow best the true orbital motion of the WR star. In 
Figure [8] we have plotted the orbital solution obtained from 
Hen A4686. Overplotted are the combined RVs obtained 
from N IV A4058, shifted by +110 kms -1 to match the mean 
velocity of the orbital solution. While the scatter of the data 
points is expectedly large (see Section [3.ip . which makes it 
impossible to use this line for a sensible orbital fit, the points 
do follow the orbital motion and, more importantly, do not 
show a significantly different RV amplitude. Even if this is 
not fully conclusive evidence, it is a very plausible indicator 
that the orbital parameters obtained from the RVs of the 
two Hell lines are reliable. Moreover, since the orbital so- 
lutions were obtained from a combined set of spectroscopic 
and polarimetric data, and the latter is not affected by line- 
profile variations, the influence of WWC excess emission is 
further attenuated. 



Since N IV A4058 does not suffers from radiative-transfer 
effects like the Hell lines, and is not blended, it also yields a 
relatively precise estimate of the true systemic velocity of the 
WR stan lMoffat fc Seggewissl 1 19791 ) found a small, negative 
shift for Galactic WN7/WN8 stars, (-15±13) kms -1 (stan- 
dard deviation). Using the mean velocity of the Niv A4058 
line as an estimate of the true systemic velocity of the binary 
system, we thus obtain W(4058) = (220 ±30) kms -1 . This 
value is just inconsistent with the report ed systemic veloc - 
ity of the LMC, v sys = (280 ± 20) kms^ 1 (|Kim et al.lll998T l. 
Either for some reason, Niv A4058 has a rather exceptional 
negative RV, or R145 indeed has a slightly negative systemic 
velocity with respect to the LMC. 

Based on our error-estimation method, the error for the 
inclination angle i was found to be 7° and 6° for the two 
lines Hell A4686 and A5411, respectively. Choosing a higher 
cut-off value than the 5% used will of course lead to a larger 
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Figure 7. Relative EWs (normalized to unity around apastron, 
i.e. cj> = 0.5) of the two emission lines Hell A4686 (filled lozenges) 
and A5411 (crosses). While Hen A5411 yields more noise because 
it is the weaker line, Hell A4686 suffers from almost twice the 
line-strength variability due to increased WWC excess emission 
around periastron (i.e. <f> = 0). Yet both lines yield consistent 
orbital parameters. 




Figure 8. Orbital solution as obtained from the He II A4686 emis- 
sion line (solid line) and RVs obtained from Niv A4058 (lozenges; 
filled symbols are from M89, empty symbols are based on the SOS 
data). While the scatter of the RVs obtained from the N IV A4058 
emission is very large, there is no indication that the motion of 
Hell A4686 is systematically different. 



error on i (and all other parameters), but the \ 2 deteriorates 
very quickly when the inclination angle varies by much more 
than 10 to 15°, and the discrepancies between the fit and 
the data become noticeably large. However, it is well known 
among polarimetrists that, due to an asymmetric error dis- 
tribution, there is a bias of the polarimetric determination 
of the inclination angle towards higher values than the true 
value of i; this bias is more pronounced, the lower the true 
inclination angle and the lower the data quality is (in terms 
of both the number of available data points and their respec- 
tive individual errors); see e.g. Aspin et al. (1981), Simmons 
et al. (1982) or, more recently, Wolinski & Dolan (1994) for 
a comprehensive discussion of this issue. Thus, our analysis 
is bound to overestimate the true value of i, which means 
that the true inclination angle is more likely lower than 39° , 



rendering the resulting masses even higher and less plausi- 
ble (see below). If, on the other hand, we were to force a 
fit with the inclination angle fixed at 90°, the resulting \ 2 
would more than double compared to the best-fit value. Un- 
less the orbital period is wrong altogether, it is thus unlikely 
that the inclination angle is that high. This is most unfortu- 
nate, since this result rules out the possibility of R145 being 
an eclipsing system, so that photometry does not seem to 
be a viable alternative to reveal core eclipses, although at- 
mospheric eclipses may still be possible. On the other hand, 
we are not aware of any existing or published photometry 
of R145, so that we have no possibility to verify the po- 
larimetric result. We will therefore adopt the average value 
of i — 38° ± 9° (quadratic error) for the rest of the paper, 
and discuss the consequences for the stellar masses further 
below. Given the great importance of an accurate determi- 
nation of the inclination angle, we strongly encourage that 
more high-quality data be collected. 



3.4 Search for the Companion in R145 

As both S08 and M89 reported, no trace of the compan- 
ion's absorption lines could be found in the individual spec- 
tra. It would be fairly obvious if the secondary were an 
emi ssion-line star as well , like e.g. in the case of NGC3603- 
Al l|Schnurr et al. 2008a). This is already an indication for 
the secondary being considerably fainter than the primary, 
and most likely an absorption-line O star. 

In a first attempt to isolate the spectrum of the com- 
panion, we shifted all spectra of R145 into the frame of refer- 
ence of the WR component, and co-added them using their 
respective S/N values as weights to obtain a high-quality, 
mean WR spectrum. This mean was then subtracted from 
each individual spectrum, thereby removing most of the 
WR star's emission lines; thus, if any O-star absorptions 
were visible, we expected to find them in the individual 
residual spectra. However, there was no obvious trace; any 
absorption- like structure was essentially lost in the noise. 

We there fore adopted a modi fied version of the method 
developed bv lDemers et al.l(|2002l ), the so-called "shift-and- 
add" method. Originally it consists in shifting all spectra 
into the frame of reference of the WR star and subtracting 
the mean (which we have done) , shifting the residual spectra 
back into the frame of reference of the O star, and co-adding 
all shifted residuals to obtain a cleaner mean spectrum of 
the O star. While this method can be used iteratively to 
obtain ever better disentangled spectra of the WR and the 
O components, respectively, its caveat is obvious: one has to 
know the orbit of the O star, otherwise one cannot shift the 
residuals into its frame of reference. 

Since we had neither the companion's orbit nor any, 
even weak, spectral signatures of absorption lines, we had 
to assume that i) the companion is indeed an O star and not 
another emission-line star (as suggested by the orbital be- 
haviour of the emission-line profiles) , and ii) the companion 
moves in perfect anti-phase (i.e. phase shift of A(/> — 180°) 
to the WR orbit. The orbit of the O star can then be re- 
constructed by mirroring the orbital motion of the WR star. 
To do so, we considered the differential motion <5i?Vij,wR of 
the WR star with respect to its systemic (zero) velocity, as 
measured at two dates i, j, and multiplied it by —q, where 
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q = Ko/Kwn = Mwr/Mo the mass ratio of the binary 
system, to obtain the corresponding 5RV[j z o- 

For instance, if the WR star changes its RV by 50 
kms -1 , and if, say, q = 2, then the O star's RV has to 
change by —100 kms -1 in the same time interval [i;j]. i.e. 
in perfect anti-phase. To reduced the influence of errors on 
the individual RVwR,(t) measures, we used the orbital solu- 
tion of the WR star as obtained above from the He II A4686 
line. We emphasize that because this is an entirely differ- 
ential approach, one does not need to know the systemic 
velocity of either the WR star or the O star (which must be 
different from that of the WR star, not least because most of 
the WR emissions display a systematic redshift; see above). 

With this tentative RVo(t) curve for the O star, the 
individual residual spectra (i.e., after WR subtraction) can 
now be shifted into the tentative frame of reference of the O 
star. Co-adding of the thus shifted residuals should make the 
O star's absorption spectrum appear with greater clarity, if 
one has the correct value for q, because the co-added aver- 
age spectrum has a higher S/N. One can then iterate the 
procedure: The mean O-star spectrum is subtracted from 
the original spectra, and the residuals are used to construct 
a "cleaner" mean WR-star spectrum, i.e. one which now is 
relatively free from any O-star absorption lines. With the 
new mean WR spectrum, the procedure to obtain the O- 
star spectrum is repeated, and so forth. 

But how can one know that a chosen q is the right 
one? If the residual spectra are shifted by exactly their cor- 
rect RVo(t) value, the co-added absorption lines will have 
maximum depths, because they are co-added while in per- 
fect superposition. If, however, a wrong mass ratio was as- 
sumed (too large or too small), then the co-added absorp- 
tion lines will be shallower, because the individual spectra 
are not shifted into the correct frame of reference in which 
they should stand still over time. Thus, when using different 
tentative values for q (or Ko, which is immediately given 
by qKwR) and measuring the depth of a given, co-added 
absorption line, the depth of the line will go through a max- 
imum when plotted versus q. The only challenge is to wisely 
choose the O-star absorption lines one measures, given that 
the WN star in R145 shows strong emission lines. 

We opted for the two lines Hen A4200 and A4542, be- 
cause WN emissions at these positions are relatively weak 
and because these are the strongest He lines in a hot, early- 
type O star, the anticipated companion. We carried out 
the "shift-and-add" exercise scanning through values of q 
so that the (intuitively more easily understandable) Ko val- 
ues range from —100 kms -1 (i.e., the O star is less massive 
than the WR star and moves in phase with it, which is of 
course unphysical, but it gives a strong constraint on the 
zero level of the measured absorption-line depth) to +600 
kms (i.e., the O star is ~ 8 times less massive than the 
WR star and moves in anti-phase with it, just as it should 
if it was its true companion), depending on the line. Thus, 
we also considered the case were the O star is actually more 
massive than the WR star, i.e. it moves with a smaller RV 
amplitude than the WR star. 

From the co-added, mean O-star spectrum, the depths 
of the absorption-lines were obtained by hand using the 
ESO-MIDAS task GET/GCUR. Continuum values on ei- 
ther side of the absorption were averaged and subtracted 
from the peak value such as to obtain positive values for the 
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Figure 9. Measured depths of the absorption lines in the co- 
added spectra versus the assumed value for Ko, shown for Hell 
A4200 (upper) and Hell A4542 (lower panel). Note the strong 
asymmetry. The parabolic fit is shown considering all data (dot- 
ted line) or only the central cusp (solid line). See text for more 
details. 



line strength, with deeper lines yielding larger values. Re- 
sults are shown in Figure [5] for the Hell A4200 and the Hell 
A4542 absorptions, respectively, where measured line depths 
versus the respective Ko value are plotted. The asymmetry 
of the curve might be caused by i) an intrinsic absorption- 
line asymmetry due to a modest P Cygni profile in the O 
star spectrum, or ii) a residual emission of the WR star 
due to imperfect subtraction. Both effects will distort the 
curve because the intensity is not the same on both sides of 
the absorption- line center. Asymmetry is particularly visi- 
ble in the data for the Hell A4200 line; however, whether 
this comes from an intrinsic asymmetry in the O-star ab- 
sorption or is related to imperfect WR subtraction, cannot 
be determined. 

In order to determine the value for Ko at which the 
absorption- line depth was maximal, a parabola was fitted 
to the data (see Figure [9]). Comparing Hell A4200 and Hell 
A4542 immediately reveals that the strong asymmetry in 
He II A4200 leads to a position of the global maximum which 
is different from the value obtained from Hell A4542, (306 ± 
69) versus (217 ± 24) kms -1 . While technically, these two 
values are consistent with each other within their respective 
errors, we re-fitted both sets of points considering only the 
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central cusp (fit shown in solid line) . We now obtained (265± 
67) and (214 ± 19) kms' 1 for Hen A4200 and Hell A4542, 
respectively. Since the error on the result for Hell A4200 is 
~ 3.5 times larger than for Hen A4542, this leads to a very 
small relative weight (~ 1/12), if one were to calculate a 
cr-weighted average. Furthermore, S08 reported that due to 
the lack of arc-comparison lines towards the blue end of their 
spectra, RV scatter was found to be larger than in the central 
region of the spectra (i.e., around Hell A4686). Therefore, we 
did not calculate a weighted mean shift but rather entirely 
relied on the value obtained from the He II A4542 absorption. 
The second iteration yielded only a marginal change, Ko = 
(212 ± 19) kms -1 . The final, resulting O-star spectrum for 
the two absorption lines is shown in Figure [lOl 

Strong artifacts are present, possibly due to imperfect 
WR subtraction, not least because the WR emission lines 
display phase-dependent excess emission due to wind-wind 
collisions (see Section |3.6[) . Moreover, upon subtraction of 
the mean WR spectrum, slight misalignments between the 
WR mean and an individual spectrum will generate pseudo 
P Cyg profiles in the residuals, whose strength and width 
sensitively depend on the amount of misalignment; at least 
partially, these artifacts are reduced by the averaging pro- 
cess. Unfortunately, however, this effectively inhibits one to 
apply this method on the Hell A4686 emission line, which 
would be much better visible than any absorption lines, 
should the secondary be an emission-line (e.g. an extreme 
Of/WN6) star as well. 

The absorption lines of the putative secondary are very 
faint; the depth of the Hell A4542 absorption line is ~0.013 
in units of relative flux, its EW = (60 ± 6) mA (error es- 
timated from results obtained for different continuum po- 
sitions), and its FWHM = (5.4 ± 0.2) A, the latter be- 
ing consistent with what is to be expected from a typical 
O star. Given the extremely high S/N obtained in the co- 
added spectra (S/N ~800), it is thus not implausible that 
these lines are indeed real. However, the low spectral reso- 
lution of our data, the frequent rebinning during the shift- 
and-add, errors in the orbital solution of the WR star, and 
noise will degrade the relatively narrow absorption lines of 
the O star, so the true depth might be somewhat larger. 
As an estimate, we adopt a degradation of the depth by 
a third, i.e. the true depth is 1.5 times larger than mea- 
sured, i.e ~0.02 relative flux. Comparing with the depth of 
the absorption lines Hell A4200 and A4542 in single, early- 
type O stars, which are typi cally 0.12 and 0.16 in contin - 
uum units, respectively (e.g.. IWalborn fc Fitzpatrick| [l990i. 
we then find that the companion is diluted by a factor of 
at least ~8; the WN6h component is thus at least ~2.2 
mag brighter than t he O star. The same exer cise can be re- 
peated for the EW. IConti fc Alschulerj (|l97lh report a typ- 
ical log£W(4542) = 2.8 ± 0.1, i.e. EW(4542) ~ 630±^° 
mA for early-type O stars, with little variation between 
the spectral types. Using the above measure, we thus obtain 
a dilution factor llijj, consistent with what was derived 
above. At ~same bolometric correction, this light ratio di- 
rectly translates into a luminosity (and hence, mass) ratio. 
Using L oc M a with a = 2 yields the mass ratio of the 
system q = M W r/M = 3.3±q 3 b . 
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Figure 10. Co-added residual spectra for Hell A4200 (upper) 
and Hen A4542 (lower panel) after being shifted into the O-star's 
frame of reference, using the best solution for Ko = 212 kms" 1 
(see text for more details). The respective positions of the He II ab- 
sorption lines are indicated. Also indicated is an unidentified emis- 
sion feature bluewards of Hen A4542, which seems to be present 
in He II A4200 as well, but much less pronounced. Whether this 
feature is real or merely an artifact of the shift-and-add method, 
is unclear. 



3.5 Masses of the Binary Components 

To determine the absolute dynamical masses of both com- 
ponents in R145, we have averaged the values of the orbital 
parameters obtained for the two Hen lines A4686 and A5411 
(cf. Table [5]). With the usual relations 



Mi sin 3 i = 1.035 x W~ 7 (K! + K 2 ) 2 K 2 P(l - e 2 ) 3/2 
and 

M 2 sin 3 i = 1.035 x 10 -7 (iti + K 2 fK x P{l - e 
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where P is in days, K in kms , masses M in solar units, and 
taking the values Kwn = (87±21) kms' 1 , e = 0.695±0.028, 
P = (158.8 ± 0.1) days, and Ko = (212 ± 19) kms" 1 from 
the Hell A4542 absorption (see above), we obtain minimum 
masses MwRsin 3 i = (116 ± 33) Mq for the primary and 
Mo sin 3 i = (48 ± 20) Mq for the secondary. Thus, the mass 
ratio is q — 2.4 ± 1.2, consistent with the value that was 
derived from the absorption- line strength (see above). 

Although the uncertainties are large, minimum masses 
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are already very high and consistent with R145's absolute 
visual magnitude. From the BAT99 catalogue, v = 12.16 
mag and b - v = 0.03. Ad opting R144's (b - v)o = -0.25 
jCrowther fc Dessardll998r i for R145, we obtain A v ~ 0.87 
mag. With a distance modulus of DM = 18.5 mag for the 
LMC and the relation v - M v = DM + A v , it follows that 
M v ~ —7.2 mag, i. e. slightly brighte r than that of WR20a, 
M v = -7.04 mag i|Rauw et al.||2007T ). 

In WR20a, which consists of two equally massive, and 
hence equally bright, components, each one of the compo- 
nents has My ~ —6.3 mag. Since the light from R145 is 
entirely dominated by the WN6h component, the bright- 
ness difference between the WN6h component in R145 and 
one of the WN6ha components in WR20 is ~0.9 mag; this 
brightness difference directly translates into a luminosity 
(and hence, mass) ratio. The flux ratio between the primary 
in R145 and one WN6ha component of WR20a is thus ~2.3; 
assuming the same bolometric correction and L oc M 2 yields 
a mass ratio of ~1.5, i.e. Mm^wNeh ~ 125Mq. A differ- 
ent exponent in the mass-luminosity relation is possible, but 
covered by the large uncertainties on the stellar mass any- 
way, as are slight changes in reddening, intrinsic colors, and 
distance modulus. 

At first glance, this result is in excellent agreement with 
the minimum mass obtained from the orbital motion. How- 
ever, now the inclination angle has to be taken into ac- 
count. The fit to the polarimetric data returns i — 38° ± 9°, 
sin 3 i = 0.23lg Hence, the true masses of the binary com- 
ponents are higher by a factor 4.3l*'g, at least 300 and 125 
Mq, respectively. Not only would such high masses be un- 
precedented, but R145 would also be severely underluminous 
for its mass. If the correct inclination angle were closer to 
60°, the masses would still be ~1.5 times higher, but then 
the resulting absolute masses would be roughly (if just) con- 
sistent with the intrinsic brightness of the system. However, 
the \ 2 °f the fit quickly deteriorates if the inclination an- 
gle is 10 to 15° larger than the optimum value (see above), 
so there seems to be no margin for a significantly higher 
inclination angle. 

Unless we accept that the polarimetry yields an entirely 
wrong inclination angle (which then would question the va- 
lidity of the orbital solution altogether) , the easiest way out 
is that the shift-and-add method is at fault, and that the 
secondary has a smaller RV amplitude, i.e. the mass (and 
hence, luminosity) ratio is smaller (A larger RV amplitude of 
the secondary would not remedy the situation, because the 
resulting system mass would be even higher, if everything 
else remained the same.) Then, however, the question arises 
why the secondary is not visible in the composite spectrum, 
let alone during the shift-and-add process. An explanation 
might be that the secondary itself is massive (luminous) 
enough to display a significant stellar wind, and resembles 
e.g. an extreme Of/WN6 star. The Hen A4686 emission of 
such a star would be weak enough to be dominated by the 
much stronger line of the WN6h component, and possibly 
even be drowned out by excess emission arising from wind- 
wind collision in the residual spectra (see Section I3.6|l ; on 
the other hand, the He II A4200 and A4542 absorptions could 
(at least partially) be filled up by wind emission. The latter 
would seriously affect the usefulness of measuring the depth 
or the EW of the absorption lines in the co-added resid- 
uals (any further challenges with this notwithstanding, see 



above), and entirely falsify the luminosity (and hence, mass) 
ratio that was derived from the absorption-line strengths, 
etc. In the following section, we will therefore use an inde- 
pendent approach to obtain an estimate on the inclination 
angle. 

3.6 Wind- Wind Interaction Effects 

As mentioned above, it is well known that in WR+O bi- 
naries, where both stars have strong winds, WWCs occur 
which generate two cone-shaped shocks separated by a con- 
tact discontinuity. This contact discontinuity wraps around 
the star with the weaker wind (i.e. that with less wind mo- 
mentum), which is usually the O star. As the shocked mat- 
ter flows with high velocities downstream along the contact 
discontinuity, it cools, giving rise to excess emissions which 
can be seen atop some optical emission lines of the WR 
star. The position of the excess bumps atop the emission 
lines depends on the orientation of the flow. The shock-cone 
acts like a beacon of a lighthouse: if it is pointing towards 
the observer, the excess emission will be blue-shifted, and 
half an orbit later, the matter will be flowing away from the 
observer, resulting in a red-shifted excess emission. Due to 
the orbital motion, the Coriolis force will slightly distort the 
cone, tilting its axis and making it curve away at large sep- 
arations, with respect to the line of sight between the WR 
and the O star. 

To study this behavior in more detail in R145, all spec- 
tra were shifted into the rest frame of the WR star. Then, a 
minimum spectrum was computed by averaging spectra ob- 
tained near apastron, i.e. around cj> ~ 0.5, when the WWC 
excess emissions are weakest. (While the excess emissions are 
not zero, as will be shown below, differences remain small 
and have no influence on the conclusions drawn here.) Due 
to the eccentricity of the orbit, phases near apastron are 
well covered, and a good-quality minimum spectrum was 
obtained. The minimum was subtracted from all spectra of 
the time series, the residuals were shifted back into the ob- 
server's frame of reference, folded into the corresponding 
phase of the binary period, and a dynamic spectrum was 
constructed where the intensities are coded in greyscales, 
white being the strongest emission. Forty phase-bins were 
used for good resolution, covering one orbital cycle centered 
on (j> — 0, at which periastron passage occurs. The result is 
shown in Figure [TT1 

Since the WWC excess emissions vary in shape and po- 
sition over the orbital cycle in a defined way, modelling their 
profile can be used to obtain an independent measure of 
the system p arameters like e.g. the orbital inclination (see 
lLuehrs|[l997l ). Unfortunately, the quality of our data is not 
sufficient to grant such an analysis, but a simplified method 
was developed by iHill et all (|2000r i that still allows one to 
obtain valuable information from the phase dependency of 
the excess emission. 

Following Hill et al.'s (2000) approach, we have mea- 
sured the full-width at half maximum (FWHM) and the 
mean RV of the excess-emission profiles, i.e. of the residuals 
after having subtracted the minimum spectrum (see above). 
These values follow the relation 

FWHM = Ci + 2v conc sin 9 sj 1 - sin 2 i cos 2 (0 - &cj>) 
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Figure 11. Ovorplot (upper panel) and greyscale (lower panel) of 
the residual spectra in the observer's frame of reference, after sub- 
traction of the minimum profile. While the He II A4686 emission 
line is highly variable, the strongest increase occurs around perias- 
tron passage (<f> = 0). The WR star passes in front at (p = 0.5 and 
it can (if just) be seen that it is at this phase that the WWC excess 
emission displays the largest redshift. Horizontal black stripes in 
the greyscale correspond to empty phase bins. 



RV 



C2 + fconc cos 9 sini cos(0 — 8(j>), 



with C\ and C2 simple constants, w CO no the stream veloc- 
ity along the shock cone, 6 the half opening angle of the 
cone, (f> the respective orbital phas^l, &4> the tilt angle of the 
cone axis due to the Coriolis force which leads to a phase 
shift between the orbital and the cone motion, and i the or- 
bital inclination angle (also see Fig. 6 of lHill et al.ll2000l ). A 
constraint is that the stream velocity in the cone must not 
exceed the wind terminal velocity (1300 kms -1 , see Section 
I3.7|l . A least-square fit to the data is shown in Figure 1121 It 
is satisfactory for the FWHM, but less so for the mean RV, 
most likely to insufficient data quality. We obtain G\ = 30 
kms" 1 , Ci = 220 kms" 1 , v coae = 810 kms" 1 , 6 = 69°, and 
5(j> — 0. The inclination angle is found to be i = 45°, in 
excellent agreement with the result obtained by polarime- 
try. Just as in the case of the polarimetry (see above), a fit 
with i — 90° is ruled out by the data; thus it appears that 
R145 is indeed seen under such a low inclination, with the 
consequences discussed above. 

As can be seen from Figure [TTl the WWC excess emis- 
sion is by far strongest at periastron passage, i.e. when 
the two stars are closest to each other, and quickly drops 
when the separation between the stars increases. The de- 



3 <j> is for a circular orbit; for an elliptical orbit, we replace <j> by 
v + u> — 90° , where v is the true anomaly. Thus, the WR star 
passes inferior conjunction when v + u = 90°. 
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Figure 12. Fit to the FWHM (upper panel) and mean RV (lower 
panel) of the WWC excess-emission profile. See text for more 
details. 



pendence of the excess emission on the sepa ration can be 
studie d in more detail, using the recipe by IStevens et al.l 
(1992). These authors have shown that the WWC zone of 
relatively wide binaries can be treated as adiabatic, even 
at periastron; therefore the amount of X-ray flux generated 
by WWC is expected to go with the inverse separation be- 
tween the stars (also see IUsovlll992l ). Since Hell A4686 is 
a recombination line (i.e. f ormed through "colli sion" of two 
particles) the formalism of IStevens et al.l ljl992T ) can be ap- 
plied to the strength of the excess emission, too. Assuming 
the WWC excess emission in Hell A4686 to be optically 
thin, we can thus expect a rise and fall of the emissitivity 
which is oc where D is the separation between the WR 

and the O star. We here use the relative separation D(v) as 
a function of true anomaly v, normalized by the minimum 
separation d m i n , i.e. 



D(v) = d(v)/d mi n 



1 + e 



1 + e cos v 



where v is the true anomaly of the orbital ellipse and 
D m i n — 1 at periastron (v = 0, phase of the elliptical or- 
bit = 0), and consequently D max = (1 + e)/(l — e) at 
apastron (v = ir, (f> = 0.5). Note that due to this normaliza- 
tion, the expression is independent of the inclination angle. 
However, while the excess emission is minimum at apastron, 
it is not generally zero. Thus, the observed EW of the Hell 
A4686 emission is at any time the sum of the EW of the 
unperturbed WR emission, Eq, and a separation-dependent 
excess flux, X /D(y), so that EW(D) = E +X /D(v), with 
Xo a constant (the maximum excess), and D(v) as defined 
above. For verification, we also fitted a steeper dependency, 
EW'(D) = E' + X' /D(v) 2 . 

From equal-weight least-square fitting, we obtained 
pairs of (Eo;Xq) = (49.2; 13.4) for the 1/D dependency, and 
(E' ;Xo) = (51.8; 11.4) for the 1/D 2 case (EWs in A). The 
data are shown in Figure [T51 with the two theoretical curves 
overplotted. While both curves certainly look equally well- 
fitting, the 1/D 2 curve yields a x 2 that is ~10% larger than 
that for the 1/D case. Since our data are rather noisy, and 
show larger scatter around periastron than around apastron, 
it is difficult to assess the true (a posteriori) measurement 
error per individual data point. Using the scatter of the data 
around apastron for this, the normalized x2 of the fit is -~2, 
but larger if the scatter around periastron is adopted. 
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Figure 13. Equivalent widths of the Hen A4686 emission as a 
function of orbital phase (i.e. separation D). Overplotted are the 
hypothetical emissivity curves; the solid line represents an emis- 
sivity rising and falling with D~ 1 , whereas the dotted line is for 
a steeper D~ 2 law. The horizontal lines indicate the unperturbed 
equivalent width, as obtained from the fit for the D~ 1 case (solid 
line) and the D~ 2 case (dotted line). See text for more details. 



3.7 The Mass-Loss Rate of the WN Star 

Independent of any issues regarding the orbital inclination 
angle, we can obtain a rough estimate of the mass-loss rate, 
Mwr, of the WN6h star in R145, assuming that the com- 
panion has a negligible mass loss. Following the prescription 
of Moffat et al. (1998; also see references therein), we use the 
total optical depth in polarization due to electron scattering 



r = 23(tt/ c Mwr(1 + cos 



) sini/[(167r) 2 m p t; 00 a]. 



where <jt is the Thomson electron-scattering cross-section, 
f c = Io/(Iwn + Io) the intensity ratio between the O star 
and the total luminosity, m p the proton mass, Voo the wind 
terminal velocity, and a = [P 2 (Mwr + Mo)] 1 '' 3 the orbital 
separation. For f c , we use the light ratio ~ 0.125 as de- 
rived above from the absorption-line strength of the O star 
companion with respect to single earl y-type O stars. From 
tailore d atmosphere analysis of R144, ICrowther fc Dessarj 
l| 19981 ) found = 1350 kms -1 . S08 report that the FWHM 
of the Hen A4686 emission line in R145 is only slightly 
smaller than in R144; assuming that both stars are sim- 
ilar enough that the FWHM of Hen A4686 can act as a 
reasonable proxy of the terminal wind velocity, we adopt for 
R145 a value of Voo = 1300 kms^ 1 . For i and r, the averaged 
values from Table [5] are taken. 

Inserting those values in the equation given above yields 
a very high mass-loss rate of Mwr = 1.0 x 10~ 4 MQyr -1 
(remarkably, for inclination angles between 40° and 90° 
this value changes barely). From tailored analysis and us- 
ing models without i r on-lin e blanketing and wind clumping, 
ICrowther fc Dessard \ 19981 ) found M = 1.5 X 10~ 4 MQyr -1 
for R144. Introducing typical clumping factors / ~ 0.1 
usually reduces the mass-loss rates derived from spectro- 
scopic means by a factor of l/y/J ~ 3, i.e. to more sensible 
M = 5 x 10~ 5 MQyr" 1 ; since R144 is the more luminous 
star of the two, one would expect an even lower mass-loss 
rate for R145. However, since polarization depends only on 



the number of scatterers (i.e. electrons), polarimetric diag- 
nostics are unaffected by wind clumping, and polarization is 
indeed observed to be that high in R145. 

Thus, assuming the models to be correct, such a high, 
continuous mass-loss rate for R145 would be sensible only 
if the wind were particularly smooth. Another way to lower 
the derived mass-loss rate would be to considerably decrease 
the orbital separation a or the terminal wind speed Voo- The 
former would require a mass-ratio closer to unity, with the 
consequences discussed above, while the latter is unlikely, 
given the close resemblance of R145 to R144. We thus con- 
clude that the primary of R145 may in fact have a very high 
mass-loss rate, subject to further confirmation. 



4 SUMMARY AND CONCLUSION 

We have combin ed prev i ously published radial velocity 
(RV) data fr o m iMoffatl (|l989l ) with data obtained by 
ISchnurr et all (|2008bl ). along with previously unpublished 
polarimetric data, for the two WN6h stars R144 and R145 
in the LMC. While the former sta r was first suspected to 
be binary by ISchnurr et all |2008bl ). the latter had already 
been identified as a binary by Moffat (1989). 

While our study could not reveal any periodicity in the 
data of R144, we have, for the first time, established the 
full set of orbital parameters for R145. The orbital period 
was found to be P — (158.8 ± 0.1) days, in contrast to the 
preliminary estimate of 25.4 days by Moffat (1989), based 
on the assumption of a circular orbit. From our analysis 
of the polarimetric data for R145, the inclination angle of 
the orbital system was found to be very low, i = 38° ± 
9°. This value was confirmed by a simplified study of the 
phase-dependent excess emission on the He II A46 86 line due 
to win d-wind collision. Following the approach of iHill et al.l 
(2000), we obtained a consistent value of i = 45°. 

By applying a modified version of the shift-and-add 
method originally developed by Demers et al. (2002), we 
were able to isolate the spectral signature of the O-star com- 
panion. We found the RV amplitude of the primary (WN) 
star to be ifwN = (87±21) kms -1 , while the RV amplitude 
of the secondary (O) star was found to be Kq = (212 ± 19) 
kms . We derived the minimum masses of the WR and the 
O component to be, respectively, Mwn sin 3 i = (116 ± 33) 
M Q and Mo sin 3 i = (48 ± 20) M . While the large mass 
ratio q = Mwr/Mo ~ 2.4 ± 1.2 is compatible with the 
strong dilution of the O-companion's absorption lines, the 
absolute masses of the two stars are problematic. If the de- 
rived, low inclination angle is correct, one obtains minimum 
absolute masses of 300 and 125 Mq for the WN and the O 
star, respectively, but such high masses are not compatible 
with the observed luminosity of R145, which puts the sys- 
tem between NGC3603-A1 (see ISchnurr et ail l2008al ) and 
WR20a i|Rauw et al.ll2007l ). Given that both NGC3603-A1 
and WR20a consist of two ^equally bright stars, while in 
R145 the WN component dominates the light, masses of 
120 and 50 Mq are more likely, but this would require an 
inclination angle close to 90°, which is not supported by 
the data. Another possibility, the mass ratio being closer to 
unity, is not supported by the results of the shift-and-add 
method. Thus, something is clearly odd in R145, and more 
and better data are required. 
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From the polarimetric data for R145, we also estimated 
a rather extreme value of the mass-loss rate for the WR com- 
ponent, of ~ 10 -4 Mqjt -1 , which is suspiciously high even 
for such a massive and hence luminous object. Analyzing 
the He n A4686 excess emission from wind- wind collisions in 
R145, we find that its strength depends on the separation d 
of the two stars as as is expected from theory of adi- 
abatic winds, but that a st eeper dependence , oc rf~ 2 , which 
was reported for WR140 bv lMarchenko et ail (|2003h . cannot 
unambiguously be ruled out. 

Given the very large uncertainties on the values derived 
from our data, our findings cannot be considered other than 
preliminary. However, given the fact that R145 seems to con- 
tain one of the most luminous and thus probably also most 
massive stars known in the Local Group, we feel that any ad- 
ditional effort is well deserved to obtain more reliable results. 
Therefore, we are currently carrying out a long-term moni- 
toring of this system, along with R144, but we strongly en- 
courage independent studies of these potential cornerstone 
objects. 
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